Equilibrium of anchored interfaces with quenched disordered growth 
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' The roughening behavior of a one-dimensional interface fluctuating under quenched disorder 

^— ■») , growth is examined while keeping an anchored boundary. The latter introduces detailed balance 

' conditions which allows for a thorough analysis of equilibrium aspects at both macroscopic and 

I microscopic scales. It is found that the interface roughens linearly with the substrate size only in 

. the vicinity of special disorder realizations. Otherwise, it remains stiff and tilted. 
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Studies of inliomogcncous interface growth have thrived in a variety of physical contexts [1] characterizing phenom- 
ena as diverse as fluid imbibition in porous media [2] and crystalline surface growth on disordered substrates [3]. It 
• is known that in those situations a small amount of disorder can severely modify the interface motion and ultimately 
^ I alter its roughening behavior [1]. This is the case of time independent but spatially random growth rates arising, for 
instance, from a quenched array of columnar defects pinning the flux lines, here playing the role of interfaces, in dirty 
high temperature superconductors [4] . Another mechanism whereby the interface character results deeply affected at 
large times is realized by anchoring conditions which suppress fluctuations at the interface boundaries. Such confined 
geometries actually occur in the unbinding of polymers from a wall [5], and also emerge as domain walls of (2+1) 
, dimensional cellular automaton models [6] , as well as in stationary nonequilibrium systems in d = 1 [7] . In this work 
we focus on the combined effect that quenched disorder and anchoring conditions brings about in the steady state 
(SS) properties of Id interfaces, at both macroscopic and microscopic scales. 
G ■ The problem is most conveniently treated in the discrete formulation of growth processes. As usual [8], here we 
' represent the latter in terms of restricted solid on solid (RSOS) configurations of heights (ho), hi, ...jhL growing 
, stochastically on a substrate of size L. Throughout the evolution, fluctuations are suppressed entirely at ho, whereas 
Q ' all other hj can increase (decrease) in two height units with substrate dependent rates ej (e^ ) [9]. Due to the RSOS 
O ,, constraints — hj\ — 1, these variations can only occur at local extrema of the interface, as illustrated in Fig. 1. 

Despite its simplicity, it will turn out that this model encompasses a disorder driven transition along with unusual 
[ roughening exponents. 

^ Following most studies [1], we concentrate on the root mean square dispersion of heights, commonly associated 

to the width W of the system. Under quenched growth this further requires the averaging of W over all disorder 
CO , realizations (denoted by the brackets below), namely 

(W^(L,<)) = ({i^ [/i,(t)-J?(t)f}'^'), (1) 

o . 

-j->; : where H{t) is the mean interface height, and the overbar indicates an ensemble average over all possible evolution 
d ' histories up to time t. On general grounds [1,10] , it can be argued that ( W{L, t) ) scales as L''f{t/L^) with a universal 
^ ' scaling function behaving as f{x) ^ x'~'l^ for x \ , whereas for a; 1 it remains constant. Consequently, for t '3> 
I the width saturates as while growing as i^^^ at much earlier stages. In what follows we content ourselves with 
' O ^ studying just the former situation i.e. the stationary regime controlled by the roughening exponent . In contrast to 
^ . the dynamic exponent z, ^ lends itself more readily for a thorough analysis under growth disorder, so hereafter we 
shall work over Eq. (1) directly in the limit t ^ oo . 

a. Particle representation - To this aim, let us first consider the conditions that the SS imposes on the above 
processes. For ease of analysis we retain the interface slopes Sj = hj — /ij-i = ±1 rather than its heights, or 
alternatively the set of occupation numbers {n} = {(1 -I- s)/2} conforming a height configuration. This enables to 
; I • exploit the well known mapping [1,8] from a Id RSOS interface to a driven particle system, i.e. hj ~ ^Ylii<j^'^i ~ 1)j 
' here specially adapted to account for both anchored and free edges as well as for quenched disorder. In particular, 
height variations at the free boundary translate here into injection and ejection of particles, as is easily visualized in Fig. 
1. In this latter representation one can readily build up a SS probability measure P({n}) satisfying detailed balance, 
regardless the specific realization of disorder considered. In equilibrium, the bulk and right boundary processes will 
demand respectively (see Fig. 1) 

e^P(ni, ... ...ni) = P( m, ... ... n^) , (2) 

e^P(ni, ...?iL-i,l) = eLP(ni, ...7ii,_i,0) , (3) 
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and therefore the SS distribution must be of the form e ' ' . Speeifieally, (2) entails a hard-eore partiele potential 

Vj = V\ +Ylii<j lii(fi/^Di whereas (3) fixes the value of the additive constant. Equivalently, by introducing the total 
particle number M = rij , we thus obtain a sample dependent grand canonical form 

P{{n}) = ^ exp(-5^y,.n, ), (4) 

3 

{Vi = 0), with a particle fugacity ^ = Y\ - ^j/^'j and a normalizing partition function Z = 11^(1 + M^"^^) • 

This provides the basic elements to compute the original average (1). In this regard and for posterior use, it can 
be easily checked that the particle densities rTj involved in the calculation of height profiles are 

(i+n^Ai)"'' (5) 

while the pair correlations needed for the analysis of height-height correlators appearing in W , result totally decoupled, 
that is. Hi rij =ninj ^ j and for all disorder realizations. However as we shall see below, there is a particular 
situation, namely /x — > 1 , for which height variables become strongly correlated. 

b. Height distribution - Before embarking in the evaluation of the double average of Eq. (1), we should determine 
firstly whether the growth of a given sample can actually produce a rough interface within the equilibrium regime (4). 
It may well happen that while the width increases with the typical substrate size, the whole system becomes actually 
smooth, as is the case of the linear profile exhibited in Fig. 1 . 

In analyzing this issue we focus attention on a more microscopic level of description such as the single height 
probability density P(/ijv) at a given location N. In turn, this requires the evaluation of an M-particle partition 
function Zm = a;„/M! constructed as 

'^m=m'' $^'e-^^i...e-^^M, (6) 

jl---3M^« 

where restricts the sums to {j^ ^ ■■■ 7^ ^ [l;^] so as to keep a constant number of particles within that 
interval. Clearly, this M-body quantity (M < ) is associated to the wanted probability at = 2M — N , since by 
construction P{2M — N) = Zm / [ YljLi (1 + /ue~^^ ) ] • To account for the hard-core interactions we build up a recursive 

relation in the particle nmnber M. Introducing the auxiliary functions g{k) = m'^ X^jLi e"*^^^' , this can be achieved 

via the following identity 

u,^=g{l)iv^_, - {M- 1) (2) , (7) 

where the latter factor is defined generically as S^_^{k) = ji^ e~'^^^ k^~^^^ '"^ ^^M-k . In particular, 

Si^j_^{2) involves constrained sums having {j ^ ji ... ^ jM-2} & [1;-^] and subtracts the unwanted terms resulting 
from (M— 1) double occupation configurations included in g(l) i-u^_-^ . Using a similar criteria, we may also infer that 
'S'm-2(2) = .9(2) ~ {M — 2) S'„_3(3) , where the last term now cancels triple occupation contributions to 5„_2(2) . 
Thereafter, iterating this reasoning down to M = 1, it can be straightforwardly verified that in terms of the j-particle 
partition functions Zj = uu-/j\ we are left with 



— 



M+1 - 



M 



g{M)+J2{-iyg{M-j)Zj, (8) 



where Zi = g{l) . Although the analytic solution of such recurrence is not reachable by standard means [11], 
nevertheless it allows for a simple numerical evaluation of the height distribution, particularly at the tails where 
statistics becomes more demanding. Despite the noisy particle potential V of individual samples it turns out that for 
M :S> 1 the distribution itself collapses towards a universal -Gaussian- form. This is evidenced by the semi-log inset of 
Fig. 1, for instance, in the case of a binary concentration p of growth rates with probability — a) (5(^ — 6) . 
Also, recursion (8) reveals that in most cases the interface is characterized by a tilted profile Hn ^ ±-/V whose local 
height fluctuations (HF) set out to be quite small, even far away from the anchored boundary ho , i.e. the interface is 
stiff, in agreement with the linear snapshot displayed in Fig. 1. On approaching certain disorder conditions however, 
Gaussian tails widen significantly (slope diminishing shown in the inset), so HF take over and the stiff regime no 
longer holds. This is illustrated by the rough profile also exhibited in Fig. 1. 
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c. Roughening criterion ~ It is therefore important to focus on these special conditions so as to ensure a meaningful 
result in the calculation of (1). Let us then consider the fluctuations of the total number of particles ( a^) = {Af^—TJ ) , 
averaged over a generic growth disorder. For equilibrium regimes, this is a representative quantity to look at, as in 
terms of the occupation densities discussed above it is immediate to check that it coincides with the height variance 
at the free edge of the interface, i.e. 4((7^) = {h\ — h^) whereas by construction it constitutes an upper bound 
for all HF. In turn, due to the RSOS constraint this variance also provides a measure of fluctuations deep inside the 
bulk. To facilitate its analysis and by virtue of the products involved in Eq. (5), it is convenient to introduce the rate 
disorder variables u,; = In (e^/ej) , in terms of which ( a^) adopts the form 



1 ^ 



cosh 



(it«.) 



(9) 



Assuming a site independent distribution of growth ratios, say with mean (u) and variance , we can thus exploit 
the central limit theorem [12] for these new variables and carry out the disorder average. Hence, it can be readily 
shown that 



>2) 



j=l j=fc+l 



^ 4s ^.f^^ j v^cosh^d) 



where aj = {n^^^_^ (1 — n^+i.^)), and k is eventually a large but finite integer {k <^ L). In particular, for a 
Gaussian w -distribution (o-^) just reduces to the integrals sum (fc = 0). The point to emphasize here is that 
whatever distribution is considered in (9), the interface fluctuations are intrinsically dominated by the above series 
of integrals. Thus, apart from a bounded quantity < otj < k/A, the eff'ect of any distribution on (cr^) will 

parallel the one caused by the Gaussian disorder. On the other hand, as long as s ^ the asymptotic analysis of (10) 
shows that the value of {u) is crucial, since for large L 

<^^^ = °' (11) 

1^ bounded, otherwise , 

which flnally provides the disorder conditions able to produce a rough interface. As was mentioned earlier, this occurs 
for fugacities such that /i ^ 1* in turn favoring strong fluctuations in the occupation numbers of Eq. (4). In that case, 
|si I < 1 and height variables become tightly correlated along the substrate because hi hj — hihj = {i — J2k=i '^1 ) > 
0, Mi < j . Otherwise, the slopes Sj —1, 1 (i.e. M /L ^ 0, 1 ), so the interface becomes tilted and hardly flexible 

since as stressed before {hf — h^) < A{a-'^) \fi . 

To illustrate the effectiveness of the above criterion for non-Gaussian cases, wc address to the binary distribution 
referred to above. Using the binomial distribution to weight all possible forms in which the growth ratios a, b may 
show up on a segment of length I , we flnd 

where 2 (pi^k = fc In a -I- (? — /c) In & . Fig. 2 shows the typical behavior of this quantity for o < 1 < 6 . In agreement with 
condition (11), there exists a critical concentration Pc = [ 1 — In a/ Inh]^^ , i.e. (u) = , for which particle fluctuations 
increase as y/L (see inset of Fig. 3), while in the thermodynamic limit diverge as \p — Pc\~^ as indicated by the inset 
of Fig. 2. In nearing Pc from below (above), this disorder-driven transition brings the interface from a tilted state with 
slope —1 (-1-1 ) to a rough phase in a continuous manner. In this respect, size effects become quite noticeable when 
passing through p^ ■ Clearly, for 6 < 1 < a the rapprochement direction is inverted, otherwise the system remains 
tilted. 

It is of interest to note here that criterion (11) implies that locally, imiform interfaces (s = ) can fluctuate stronger 
than disordered ones. In fact, for = Ci all interface conflgurations are equally likely so the single height probability 
density derived from Eq. (6) reduces to P{h) ~ ^ exp(— |jj) . So, in the pure substrate oc L , rather than diverging 
as in (11). What remains to be seen is whether different scaling behaviors can also emerge on a more macroscopic 
description, such as the width of Eq. (1) and towards which we flnally turn. 
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d. Width behavior - Let us then consider this global quantity in the slope representation discussed throughout. 
After some algebraic steps, it can be readily demonstrated that the width of a given sample is expressed as 

^'^^^ = + 5^ (i + 1 - i) (i - 1) sisy . (13) 

i<j 

In particular, for equilibrium regimes there is a key simplification in the average oiW, as SiSj — 'SiSj over all disorder 
realizations. In passing, it is worth recalling that in a uniform unanchored system W ~ y'L/12 [8], thus for = 
where s « = , the anchored boundary just modifies a mere amplitude. 

For binary rates, a closed expression of (siSj) can be easily found upon recurring once more to the binomial 
distribution. According to the occupations (5), the contributions of SiSj to its own quenched average are simply 
obtained by weighting the joint number of possibilities in which the growth ratios a, b appear along the intervals 
[i,j — 1] and [j,L] . Specifically, this yields 

k=0 k'=0 ^ / \ / 

X (1 tanh^j-fc tanh6',,fc+fc/ , (14) 

where k {k' ) is the number of a-occurrences within [j, L] {[i,j—l]), while the ^-arguments are generically defined as 
^Oj^k = klna— {j + k — L — l)ln6. The size dependence of the resulting quenched width is plotted in Fig. 3 (solid 
lines), using the critical concentrations Pc required by criterion (11). Clearly, the results support a power law growth 
of ( W) consistent with a rather unusual roughening exponent ~ 1 extended over more than three decades. Thus, 
we see that in sharp contrast to what occurs at the microscopic scale where HF of disordered interfaces are weaker 
((ct^) cx ^/L) than uniform ones (cr^ L), macroscopically the former can grow rougher than the latter. 

Notice however that on approaching weakly disordered regimes size effects become quite pronounced. Even for 
large substrate lengths (W) grows with the ( = 1/2 exponent characteristic of homogeneous systems [8,13] though 
asymptotically the linear scaling is recovered. This suggests a rather abrupt change separating the uniform case from 
the more general disordered situation, namely: at the static level of the C -exponents there is a discontinuous scaling 
regime at Pc = 0+, 1~ where e ~ e' for most sites. In fact, the numerical estimation of (W) via Eqs. (5) and (13) 
over a Gaussian growth distribution satisfying ( In (e'/e) ) = , as required by our criterion, yields the same value of 
C while also implying a discontinuous behavior near the homogeneous case, when s — > 0"*" . This is indicated by the 
dashed lines of Fig. 3 after averaging over 10'* samples. Whether such discontinuous feature corresponds to a general 
aspect of weakly disordered roughening processes in = 1 seems difficult to elucidate rigorously [14]. 

To summarize, we have constructed a recursive relation [Eq. (8) ] which allowed us to examine disordered interfaces 
at microscopic scales. For most situations this yielded a stiff picture given the small HF obtained, though in the 
vicinity of special disorder realizations HF take over and wipe out the interface stiffness. This led us to propose a 
roughening criterion by looking at the fluctuations of the total particle number along with their divergence conditions 
[Eq. (11) ]. Within such regimes we analyzed the scaling behavior of quenched widths using both Gaussian and non 
Gaussian disorders. The corresponding results support a linear growth of ( W) with the substrate size, i.e. C ~ 1 
(Fig. 3), which possibly reflects a tendency of anchored interfaces to crumple on large scales [1,14]. 

Finally, though the non-equilibrium dynamics is out of the scope of this work, one could expect that criterion (11) 
will introduce somehow a crossover between ergodic and non-ergodic evolutions. It is then natural to ask whether the 
fundamental scaling between length and time - embodied in the z exponent ofW- would be aflFected by our roughening 
condition. So far, this remains an open issue. 
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FIG. 1. RSOS interface anchored at ho, thought of as 
an asymmetric exclusion process in which height slopes 
— hj — hj-i are associated to particle occupations 
nj = i(l -|- Sj). At the free edge hL variates in 2 (—2) with 
rate ei, (e^) by particle deposition (evaporation), whereas for 
1 < j < L height changes correspond to left (right) particle 
hoppings with rates ej (e^). Lower panel: typical snapshot 
for L = 10'^ after 10 ' steps per height using a binary disorder 
with e'/e = 0.8, 1.2 under condition (11), i.e. for p = pcin the 
text. Otherwise, configurations become tilted (dotted line). 
Inset: Gaussian distribution of heights for j = 100, at and 
slightly below Pc (upper and lower lines). 
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FIG. 2. Particle fluctuations arising from a binary growth 
distribution. Data for b = 1.5 while varying the concen- 
tration of a = 0.5 . Curves from top to bottom denote in 
turn results for L = 10^,10* and 10^. The inset exhibits 
the same algebraic divergence (dashed slope) either nearing 
Pc = [ 1 — In a/ In 6 ] from below or above (dotted and solid 
lines closely following each other; L = 10^). 




FIG. 3. Quenched interface width scaling linearly with the 
substrate size. Averages were taken over Gaussian and bi- 
nary distributions (dashed and solid linos respectively), using 
condition (11). Curves in descending order stand in turn for 
s = 1, Pc = OA, s = 0.01 and pc = 0.004 (a = 0.5); the 
lowermost slope is 1/2. The inset displays fluctuations (12) 
scaling as L^^^ (dashed slope) for pc = 0.4. 
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